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Field strength correlators are semi-classically evaluated in the dilute gas mo- 
del of non-Abelian sources (instantons) and compared with lattice data for 
QCD at zero temperature. We show that one of the Euclidean invariant, 
tensorial structures vanishes for configurations being purely selfdual or anti- 
selfdual. We compute the invariant functions contributing to the correlators 
within the two lowest orders in an instanton density expansion. Fitting in- 
stanton size and density for quenched and full QCD, we obtain a reasonable 
description. 
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I. INTRODUCTION 



In the last years a systematic description of non-perturbative effects in QCD has been 
given in terms of the gluon field-strength correlators They are of immediate phe- 

nomenological importance in the model of stochastic confinement of color charge, giving 
a detailed description of the level splitting of heavy and light QQ bound states, and in 
the description of high-energy hadron and quark- (anti) quark scattering |Q. By now, there 
exist numerical results from lattice simulations concerning the fundamental field strength 
correlators for pure gauge theory with gauge groups SU{2) |]5|-|7| and SU{?>) P|-pT| over 
physical distances ranging up to 0(1) fm. The correlators have also been calculated near 
to the deconfinement transition in pure SU{?>) gauge theory 0. Very recently, this study 
has been extended to full QCD with four flavours of dynamical staggered quarks JlOl- In 
this latter work, an interesting pattern of the quark mass dependence both of the gluon con- 
densate and of the gluonic correlation length has emerged. The correlation length smoothly 
changes from the lightest quark mass to the quenched results. The extraction of the gluon 
condensate itself is afflicted with some reservations about the renormalon ambiguity in the 
determination of the Wilson coefficients of the operator product expansion. In Ref. [jlO[ this 
problem is finally neglected. As demonstrated there, the emerging mass dependence can 
be consistently described by a low energy theorem which relates it to the zero-mass case 
through the quark condensate. The real problem hard to understand is the much stronger 
gluon condensate for pure Yang Mills theory. 

In the direct computations of the gluonic field strength correlators on the lattice the non- 
perturbative behaviour has been extracted with a cooling procedure which serves the purpose 
to erase short-range perturbative fluctuations. Then at various distances d (in lattice step 
units) gauge field structures are expected to become visible in the correlator if an appropriate 
number of cooling iterations {n oc d? due to the diffusive nature of cooling) is applied. In 
view of these results it has become very likely that the field correlators have a semiclassical 
origin and might provide information about the importance and the physical parameters of 
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those classical configurations forming the basis of the semiclassical approximation. 

Instantons - solutions of the Euclidean Yang Mills field equations - are well-known ex- 
amples of semiclassical configurations, which exist both in the continuum and on the lattice. 
Their contribution to the path integral quantization plays an important role in explaining 
chiral symmetry breaking and many phenomenological facts (see e.g. [0-|TB[). 



It seems quite in time to compute the field strength correlation functions within a semi- 
classical model based on instantons and to point out to what extent the latter are able 
to describe the results of the lattice measurements mentioned above. In principle, such a 
comparison should enable us to quantify the importance of the semiclassical modes forming 
the QCD vacuum or finite-temperature Euclidean fields, respectively, and to provide their 
relevant physical parameters. 

Performing the analysis of field strength correlators with the simplest instanton solutions 
we will restrict ourselves to the SU{2) case, being aware that the embedding into a larger 
SU (Nc) gauge group can be easily realized. Therefore, the parameters of instanton densities 
and sizes we are going to determine from a comparison with lattice data will always refer to 
the physical SU (3) case. 

Earlier studies of field strength correlators due to instantons can be found in |]TP|,^ 



There the results for the single instanton approximation were obtained in terms of two- 
dimensional integrals. Here we pay particular attention to the path dependence of the phase 
factors. This can only be done numerically which was beyond the scope of these papers. 



Ref. pO| points out the insufficiency of the single instanton approximation. We study here 
also the next order terms in a density expansion. With the numerically evaluated invariant 
functions, we estimate instanton density and size. 



II. THE FIELD STRENGTH CORRELATOR 
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A. General Statements 



The gauge invariant two-point correlators of the non-Abehan field strength are defined 

as 

Dfip,ua{xi -X2) = {0\Tr{Gf,p{xi)S{xi,X2)Gi,a{x2)S\xi,X2)}\0) , (1) 

where G^p = T"'G'^^p is the field strength tensor and S{xi,X2) is the Schwinger-line phase 
operator, i.e. the parallel transporter necessary to join the field-strength operators at points 
xi, X2 in order to respect gauge invariance. denotes the generators of the gauge group 
SU (Nc). The most general form of the correlator compatible with Euclidean 0(4) invariance 
at zero temperature is, in the notation of Ref. [0], 

Dpp,uaix) = {6pjp^ - 5p^5py) {d{x^) + Di{x^fj + e^p^^ D2{x^) + 

(2) 

,dDi{x^) 

~\~ [X nXi/d pu X uX(jO py ~t~ X pX(jd pi/ X pXyO u 



„2\ n /^2\ n^/^2\ 



with X = Xi — X2 and D{x ) , Di{x ) , D2{x ) representing invariant functions. The 
invariant i52-term has been added for later use. It is only relevant for cases in which GP 
symmetry is violated or a restriction to sectors with topological charge Qt ^ is possible. 
It has been shown that Di does not contribute to the area law of Wilson loops . In 
the perturbative regime both invariant functions D and Di behave like . Only 

Di receives a contribution from one-gluon exchange |l|J^. On very general grounds the 
perturbative contribution to D (which appears at one loop and higher orders) was recently 
shown to be cancelled by higher correlator contributions Here, we shall not discuss the 



perturbative contributions in more detail. Instead we will concentrate on the contribution 
from semiclassical configurations. 

We consider the correlator of Eq. (^ as an average of the expression in brackets over 
gauge field configurations with a weight equal to exp(— S^;) {Se is the Euclidean action of 
the gauge field) times the determinant of the fermionic Dirac operator. We are going to 
compare with lattice simulation data which are taken partly in full QCD, partly within 



the quenched approximation when the fermionic backreaction on the weight is neglected. 
In our estimates only density and size of instantons play a role which are - in the case of 
full QCD - influenced by the effect of the fermionic determinant. The instanton solution 
of the Euclidean field equations is taken with its algebraic form in the evaluation of the 
field strength and of the Schwinger lines. The resulting expressions are then contributing in 
leading order to the functional average. 

At this place let us make a general statement concerning the leading order contribution of 
pure (anti-) self dual gauge field configurations to the field strength correlators. This property 
of the one-instanton contribution was mentioned before in Ref. ||2^ . 

Lemma. Pure (anti-) self dual field configurations can contribute only to D and 1^2, but do 
not contribute to Di. 

In order to prove this we assume that the gauge potential A satisfes the (anti-) self duality 
condition on the field strength 

Then the duality transformation ^Sa-yiMp^Spsi^a , acting on both field strengths in Eq. (^, 
provides 

(TrGa^SGf^sS^) = {D{x'^) + Di{x^) + 1 {^o.p^'yS - ^aS^p-!) + £^7/35 D2{x^) 

(4) 

dDi{x^) 



- {XaXf^S^s - XaXsS/3^ + X^Xs^a^ - X^X^das) 



dx 
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The sign of the third term of the right hand side of Eq. @) is just opposite to the sign 
of the corresponding term in Eq. (^. The left-hand sides of (0) and (H) are equal due to 
(anti-)selfduality. Thus, we get Di = const . From the requirement that the correlators 
should vanish as ^ oo we arrive with Di = , q.e.d. 

Corollary. Consequently, a nonzero Di may appear only due to perturbative fluctuations 
and in higher order of a density expansion for distinct semiclassical contributions (due to 
the non-linear interference of pairs of (anti-)instantons). 



B. Lattice Results 



The field-strength correlators (|^) have been estimated from lattice two-point functions 
for pure SU{2) and SU{3) gauge theory at T = and T 7^ as well as for full QCD 
PJ5|J^^T(1|. After applying the cooling method the non-perturbative contributions to the 
correlators have been shown to fall off exponentially at distances between 0.1 and 1 fm. Fits 
to the cooled Monte Carlo data still exhibit a remnant of the perturbative tail ~ x^'^ . 

In the following we want to consider the zero temperature case which, in the lattice 
measurements, is described by the following two correlators: 



dDi{x^) 



dx'^ 

(5) 

D^{x^) = D{x^) + Di{x^) . 

We will refer to the best fits to pure SU{3) gauge theory data obtained for several bare 
lattice couplings in M: 



D{x ) = Aexp{-\x\/XA) + ^exp(-|a;|/A« 

X 



b 

Di{x^) = B ex-p{-\x\/\B) + exp(-|x|/A6 

X 



(6) 



with 



A o B 



,,.-3.3x10^, -x^0.7xl0«, a^O.69, 6^0.46 
Af Af 



In physical units, obtained from string tension data, this amounts to A^i ~ 0.22 fm, ~ 0.43 

fm. The fit was achieved with a very reasonable /N^.f. — 1.7 . 
There is another solution reported in ^ with 

-x -2.7x10^, TT~0. , a ^0.4, 6^0.3, 

= Ab ~ -r-T^ , Aa = Af, ^ 00 . (8) 

1\L -Loo 
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and a similar value. 

The results show the invariant function Di(x^) to be by a factor 0(5) smaller than 
D{x^) . This indicates that D{x^) is mostly of non-perturbative nature. For the per- 
turbative contributions alone the opposite relation Di > D has been proven to be valid 

Within full QCD the same field strength correlators have been measured for Nf = 4 
flavours of staggered fermions with the Hybrid Monte Carlo (HMC) method |]lOl- Results 
were obtained for two quark masses = 0.01, 0.02 (in lattice units). Fits with the same 
ansatz (^ and (^ have provided, for the case rrig = 0.01, the following parameters 

A B 

-^^1.74x10^°, -^^0.2x10^°, a^O.71, 6^0.45, 

with a x^/^d.f. — 0.5 . Determining the scale parameter Ap for full QCD from the estimate 
of the p-meson mass one has Xa ~ 0.34 fm, which is somewhat larger than the quenched 
value quoted above. For the other quark mass value ~ 0.29 fm was obtained. We 
conclude that qualitatively the functions D and Di behave similar in all cases considered. 

In a recent paper [|ll|] the field strength correlators have been investigated in the presence 
of a static qq pair for quenched QCD (see also 0). The exponential decay of the form 
factors D and Di has been seen in this case even without cooling. The correlation lengths 
estimated are compatible with the results mentioned before. 



C. Semiclassical Approximation 

Next we want to work out the lowest order semiclassical approximation for the field- 
strength correlator by separating the Gaussian integral over quantum fluctuations on top of 
a single classical background field from the gauge invariant product of field strengths simply 
evaluated for that classical field. What then remains are zero-mode integrations over the 
appropriate set of collective coordinates F characterizing the classical field with a density 
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function M{T), 

D,p,U^) = ^J d^M{^)■Tr{G^,{x^■,^)S{x,,X2■,T)G,Ax2■,T)S\x2■,T)}, (10) 

where G^u{x; T) is the field strength tensor corresponding to configurations A^{x, T). To be 
more specific, we imagine a model of the vacuum state that is semi-classically represented 



by superpositions of instantons and anti-instantons [|l3 

N N 

Af,{x, r) = ^ a^(x; 7^) + J2 ^i^i^'^ 7i) • (H) 
i=i j=i 

The 7i ( 7j ) denote the collective coordinates of the z-th instanton (j-th anti-instanton), 
which include the position Zi, the group space orientation Ui, and the size pi. The integration 
measure in Eq. ([1^) is then expressed by 

N N 

dT = Y\_ ^li n ^Ij ) "^Ti = d'^ZidoJidpi , d'^j = d'^Zjdujjdpj . 
i=\ j=i 

For practical use we will consider here only (anti-)instantons of fixed size. This corresponds 
to the instanton liquid model invented in with a delta like size distribution. In principle, 
a more realistic p-distribution with a selfconsistent exponential infra-red cutoff (allowing to 
satisfy low-energy theorems) can be obtained from the assumption that (anti-)instantons 
repel each other at short distances [0,0. 

If the instanton liquid or gas is sufficiently dilute we can approximate the functional 
integral by an expansion in powers of the (anti-)instanton densities 714 = N/V ( 724 = N/V ). 

Strictly speaking, the superposition ansatz ([TT|) makes sense as an approximate saddle 
point of the action only if the vector potentials a^, cifj, decrease fast enough. This happens 
when the singular gauge expression is used for the (anti-) instanton solutions a^, d^, instead 
of the regular gauge form W^. The existence of a systematic expansion in higher order 



contributions to the measure A4{T) has been proven in Ref. It will be pointed out in 
the next section that, as long as one remains within the single instanton approximation, the 
actual choice between an instanton in the regular or singular gauge does not matter if the 
field strength correlator is evaluated on this background with a straight line Schwinger phase 
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factor inserted. In the calculation of interference terms containing a^(a:;7) and a^(x;7), 
however, the choice of the singular gauge form of solution is essential, and will be used in 
section 0. 

III. FIELD STRENGTH CORRELATOR IN THE ONE-INSTANTON 

APPROXIMATION 

Let us compute next the explicit single-instanton contribution to the field strength cor- 
relators (|l|) . The leading term is given by the sum of single instanton (J) and anti-instanton 
(7) contributions 



flp,U(T 



Dip,uAxi, X2) = D^p ,^^ + D^^^ 

= n4 y" d^z Tr {G^,p{xi] ^)S{xi, X2\ l)G^„{x2\ -i)S\xx, X2\ 7)) + (12) 
+ (n4,7 ^ ^4,7) • 

The integration over the (global) group orientation is trivial in this case and has been 
omitted. We emphasize, that the Schwinger line phase operator is a path dependent matrix 
in the fundamental representation 

S{xx,X2\z^=P exp y dtXp(t)ap{x{t); z)^ , (13) 

where the vector potential = T"a°, in leading order of an expansion in the density, 
belongs to the single instanton source localized at z. Together with its adjoint, , S 
takes care of the parallel transport of the field strength tensor G^p from one point of the 
measurement to the other. 

We begin with the instanton solution in the so-called regular gauge. Its SU(2) vector 
potential is expressed with y = x — z a.s follows 

a^(x;2;) = 2r7„^^-2q^ , a = 1,2,3, (14) 
G%ix; z) = -4 r]a^^ ^ ^2^2 ' = 1, 2, 3, 4 . (15) 
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The rotational degrees of freedom in the group space have been omitted. The t'Hooft tensor 



rja/j^iy (and f]a^,y for the anti-instanton) is defined in [|T^]. The Schwinger hne phase factor 
depends on a particular path between the points of measurement xi and X2- For a straight 
line path between xi and X2, the Schwinger line can be written explicitely 

S{xi,X2;z) = P exp (iT''r]af,u f dt ^^'^ j ^ ^^g^ 

V -^0 P +y J 

where are the Pauli matrices and y = Xi + {x2 — Xi) t — z denotes a point running from xi to 

X2 relative to the instanton center z. The Schwinger line can be conveniently parametrized 

with zi = xi — z and Z2 = X2 — z as follows 

S'(xi,X2;2) = exp Qr^n^e^ . (17) 



where 



Q = 2(zlzl-{z,Z2f)'''^ , ^ = arctan(x.)-arctan(xi) 

" ^ {{z,-Z2f p^ + zlzl-{z,Z2ff ^ ' 



and 



Zi{z2 - Zl) 
{{Zi-Z2f p^ + zlzl-{ZiZ2n 



y. = "^^^2 ~ 1 = 12 (19) 



Here denotes the components of a unit vector in isotopic space, 



(2:2 ^2 _ (^^^2)2) 



In what follows we shall compare the straigth line case ([T7|) - ( |20D with a special case when 
the Schwinger (path dependent) phase factor becomes trivial as S{xi,X2; z) = 1 for all points 
Xi and X2. This happens if the path connects the points Xi and X2 via the instanton center 
z along the two radial rays. In this case it is very simple to perform the integration over z 
in ([T2|). The result can be expressed by the following function of x/p = \xi — X2\/p 

( p) = / "^''WTTnWTTf • ^^^^ 



The result of the integration in (|T2D for an instanton is 
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(22) 



where 



\p 2 x^\c 



and c = + 4p^x^, which coincides with the result obtained in 
behaviour of this function is 



hm /r ( — I — > — , 



hm Ir ( — 



In - . 



I (23) 

The asymptotic 

(24) 



p J x-^QO J X p 

The correlator due to an anti-instanton is obtained by the replacement i] ^ fj which results 



m 



Upva } 



8 n^lr 



X 



(25) 



In the instanton-anti-instanton dilute gas the observed correlator in the lowest order of the 
density expansion is the sum of (^) and (P^f). Under the standard assumption that the 
instanton and anti-instanton scale sizes are equal we obtain for the invariant functions 



Dix' 



X 



(^4 + n4)/r 1-1 , 



DAx^ 



0, D^ix' 



x 



{n^-n^)Ir ( - j . (26) 



The vanishing of the one-instanton contribution to Di illustrates the general theorem 
formulated in section [11 A| . Only in the case that the densities n4 and are different, 
the e tensor structure of the field strength correlator is different from zero. In lattice 
computations one could select gauge field configurations according to their net topological 
charge Qt ^ ^ ■ Such an ensemble would allow to extract D2{x) . However, this is not 
the standard case. Therefore, we shall put = 77,4 in the following. 

If the Schwinger line is chosen to connect the points Xi and X2 along a straight path some 
corrections to formula (p3D appear. One has to replace in ( p3D Ir{x/p) by Iph{x/p) where 
V(x/p) = Wp)+^^h{x/p) 



d'^z cos 9 



P 



[zi+p^izi+pr ' 



(27) 
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with G given in ([T8|). The integral /i cannot be evaluated analytically. Let us estimate the 
degree to what the field strength correlators depend on the choice of path in the Schwinger 
line. For the straight-line path between the points xi, X2 the result can be rewritten from 
(0), (0) as 

where B is given in (^) and the average is with respect to the measure in (|2T|). An estimate 

for G can be given for x = xi — X2, with the instanton center located at distance R (with 

X + X 2 \x\ R (G^) 

R ^ \x\) from the midpoint ^ 2 — ~- © — j^j _|_ — i- Thus the difference 1 — (cosG) ~ ^ 2 — 

^JW^+fF^"^ — 0-2; where the typical instanton gas parameters R 1 fm, a; ~ p ~ 0.3 fm 

have been inserted. The result of the numerical integration for the straight line is compared 

with the expression ( p3D for the radial path gauge transporter in Fig. 1. The correlator with 

the straight Schwinger line is roughly 20 % smaller than for the 5* = 1 case in accordance 

to the estimate, and correspondingly the half width is smaller by roughly one third. In 

the estimation of the (anti-) instanton density and size below, the numerically calculated 

correlator with this minimal, straight path Schwinger line should be used. 

Considering the so-called singular gauge instanton solution 

2 

al{x-z) = 2nat.u^rT^r-^ ^ a = 1,2, 3, (28) 

y [y + p ) 

f-^a { . \ ^ P \^ - 1- yi^yfj. - yKyu 1 ^ o a (oq\ 

(r + p ) 12 r r J 

we notice that it is related to the regular gauge (anti-) instanton by a singular gauge transfor- 
mation. With the straight Schwinger line inserted, the correlator is identical to the correlator 
numerically evaluated for the regular gauge instanton. Just in order to see how the choice 
of the Schwinger line infiuences the field strength correlator, we can consider also a path for 
which the Schwinger line non-Abelian phase factor is equal to S = 1. This is the path leaving 
Xi along a radial ray (starting from the instanton center) towards infinity and approaching 
X2 from infinity along another radial ray, with an arc at infinity in between. 

The result of the integration over instanton position can be expressed as a function of 
x/p by the replacement I Is which is defined as 
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There exists a closed expression of the integral (see the lower curve in Fig. 1) 



+ 



1 + ^-31 + ^1 + 3- 



J \ P^ ) \ P^ ) x^ + 

The asymptotic behaviour of this function is 



(p2 + x2)v/c + x^ + 3pV 



hm /J - ^ — , lim J J - ^ ^ . (32) 

Note that in the next to leading order (i.e. the approximation quadratic in the density) only 
instantons written in the singular gauge can be employed to form superpositions according 
to Eq. (0). 

The long-distance asymptotics of the field strength correlator evaluated non- 
perturbatively with single instanton contributions resembles a perturbative contribution 
and does not show an exponential fall-off. This indicates that only a strongly interacting 
instanton gas or liquid might mimic the correct infrared behaviour of the theory. Neverthe- 
less, it is reasonable to ask, whether single instanton contributions can describe reliably the 
behaviour of the field correlators at intermediate distances, where lattice data are available. 
Independently of its concrete form the instanton contributions to the correlation functions 
should be compared only with the non-perturbative part of the lattice result, i.e. with the 
pure exponential terms in Eq. (^. 

As already mentioned the lattice data for Di are definitely smaller than those for D. 
Due to the ( ant i-) self duality, the tensor structure related to Di even strictly vanishes in 
the leading order of the density expansion in a dilute instanton gas picture. Therefore, the 
leading instanton result points qualitatively into the right direction. 

The lattice measurements of the field strength correlators are obtained with straight 
Schwinger line gauge transporters inserted between the points Xi and X2. Comparing with 
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the lattice data of Di Giacomo et al. [§] we can roughly estimate the instanton gas parame- 
ters. Fitting the instanton results for D according to Eq. ( p6D (with 1.^ , Iph , Is related 
to different choices of the Schwinger line path) to the first term of D in i.e. 

within the range 1 < x/p < 5, we obtain for p and n\ = 71^ + n^, the radius and the total 
density of pseudo-particles, the following results: 



quenched QCD 


I, 


Iph 






0.78 


1.35 


2.12 


n\ /fm~^ 


6.19 


4.03 


2.29 


p /GeV"^ 


0.87 


1.51 


2.36 


n\ /GeV^ 


9.39-10-3 


6.11-10-3 


3.47-10-3 


t 4 

nlp^ 


5.37-10-3 


3.14-10-2 


1.08-10-1 




7.51-10-2 


4.89-10-2 


2.77-10-2 



where G2 denotes the gluon condensate and n\p'^ the packing fraction. 

Strictly speaking, only the column denoted by Iph with this function evaluated in the 
instanton field according to the straight Schwinger line prescription should be compared with 
the result of measurements obtained on the lattice. We provide also the other estimates in 
order to get a feeling about the sensitivity with respect to the choice of the Schwinger line 
path. 

The difference to the value of the gluon condensate G2 = 0.14(8) GeV^ extracted for the 
quenched theory in Ref. can be explained by the fact that the latter result relies on the 
validity of the exponential part of the fit down to x = 0, while the single-instanton correlator 
itself is fiat for x 0. The instanton size, practically identified with the correlation length, is 
as usually adopted, p ^ | fm within the parametrizations leading to Iph and Is- The packing 
fraction n\ p^ almost coincides with the value for the instanton liquid. Yet, the density 
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estimated is bigger compared with the density of 1 fm~^ usually adopted in phenomenological 
applications. 

However, the instanton liquid phenomenology makes references only to full QCD with 
realistic quark masses. Therefore, we give here a fit based on the instanton gas formula to 
the full QCD lattice data, too, described by expression (^) and parameters @. We obtain 
from a fit in the range 1 < x/p < 5 the following parameter values: 



full QCD 


Ir 


Iph 


Is 




0.78 


1.28 


2. 


n\ /fm~^ 


0.72 


0.55 


0.64 


p /GeV"^ 


1.34 


2.21 


3.45 


n\ /GeV^ 


1.09-10-3 


8.34-10-^ 


9.7-10-4 




3.56-10-3 


1.97-10-2 


1.37-10-1 



Through an estimate based on the Banks-Casher formula (see e.g. the recent review by 
T. Schafer and E. V. Shuryak ||T7| ) 

{qq) = - — J^n4, (34) 
np V 2 

we obtain for the quark condensate: 



full QCD 


Ir 


Iph 


Is 


experiment 


{qq)/ GeV3 


-1.66-10-2 


-8.82-10-3 


-6.1-10-3 


-1.06-10-2 1.66-10-2 



We conclude that the nonperturbative tensor structure D{x^) in the field strength cor- 
relator at zero temperature can be roughly described by a semiclassical picture based on 
instanton-like non-perturbative field configurations. The extracted parameters are in the ex- 
pected ballpark when the fit is applied to the lattice measurements with dynamical fermions. 
This lends support to the conjecture that ( ant i-) self dual configurations are dominantly con- 
tributing to the correlator and the gluon condensate as known from real QCD. 
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The gluon condensate in the quenched theory obtained in Ref. 0, however, is almost 
one order of magnitude bigger than that for full QCD with light Kogut-Susskind quarks. 
With our instanton shape for the correlation function we can reduce this to a roughly half 
as big estimate. Still this results in an unexpectedly high instanton density. Due to the 
larger correlation length with dynamical light quarks, according to our fits, the packing 
fraction is bigger only by some 50 % in pure Yang Mills theory compared to QCD. In the 
following section we will study the influence on estimated instanton density and radius for 
the quenched case when corrections of second order in the density to the field strength 
correlator are taken into account. 



IV. SECOND ORDER DENSITY CONTRIBUTION 

In this section we shall present some estimates of the next order term in a density 
expansion. We have to consider the field strength for a superposition of solutions a and 6, 
where both a and h can represent an instanton or anti-instanton, 

G^,p{a, h) = G^p{a) + Gpp{b) + AGpp{a, h) , 

(35) 

/\Gpp{a,h) = -i{[ap,bp] + [bp,ap]} . 

We neglect for the purpose of this estimate the Schwinger phase factors which are known 
to give a 20 % effect in first order in the density approximation. The interference between 
classical solutions contributes terms of second order in the density, and the resulting field 
strength is neither selfdual nor anti-selfdual. Therefore we expect Di to receive the lead- 
ing contributions in this order. The functional weight (|T0|) specified for superpositions of 
Euclidean solutions is approximated by an uncorrelated ansatz in terms of the single- (anti- 
) instanton density. Therefore the second order contribution to the correlator has the form 

a,b=II 

(36) 

Tr {G^p{a{xi, 71), 6(xi, 72)) G^^{a{x2, 71), b{x2, 12))} ■ 
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When considering superpositions, one has to take the single (anti-)instanton field configu- 
rations a and b in the singular gauge. With the notation y — x — z their vector potential 
is expressed as follows 

2 p2 

where cUaa' denotes the relative color orientation of the pair (replace 77 —> 77 for anti- 
instantons) . 

As was mentioned above, within the approximation linear in the correlator Di is 
strongly zero. However, the interference between (anti-)instantons generates non-trivial 
contributions to it. In general, one can easily derive the contributions to the field strength 
correlators by averaging over the color orientation of the (anti-)instanton a relative to 
the (anti-)instanton b . Only the AG^p(a,6) AGi,a{a,b) contribution is non-zero after this 
average has been taken. Introducing the tensorial decomposition of the integral 

/ A vAv + ^Uiv. P)f(y + P) = s^M^. P) + vM.. P) (38) 

in terms of invariant functions Ji and J2, after some minor algebraic manipulations one can 
obtain the following system of equations for the second order terms, proportional to the 
number density of different pairs in the gas: 



Q^2 ^~ Jy2_^ 2 (3'^l(^>P)+^ J2{x,p))J2{x,p). 

The integrals Ji and J2 can be computed analytically, 

M^,p)-^^{4^ + 2^1n^-4(l + ^yin(l + ^) 

p'{V~c-x') \ 



(39) 





^3/2 




P' 




P' i 


3 


X^ \ 






7 


[7~ 



(40) 



In 



(p2 + x2)^ + a:^ + 3pV 



2 fi 2 / 2\2/2 \ / 2^ 



^2(^,p) = -^<!-8^ + 2-^ln-^-4 (^1 + -^J (^-^-2jln(^l , 



(41) 



p^ (p2-Fa;2)Vc + a;'^ + 3pV 
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where the notation c = + Ap^x^ was used again. The asymptotic behaviour of these 
functions is 

hm Ji(x, p) vrV^ , hm Ji(x, p) ^2n^^ . (42) 

27r2 n4 

hm J2(x,p) ^ — In(xVp^) , hm J2{x,p) -4:71^^ ■ (43) 

x^O 3 x^oo X 



The solution of the system of equations ( PUj) can be written as 

D(^)(x) = 9vr^ ^ ^^^i±^il^/(^)(x/p) 



(44) 



where the functions if^ have been determined numericahy. These functions turn out to 
be non-negative for ah x. Their values at zero distance are /(2)(0) = 1.28 and ^(0) = 0.28. 
Taking the first and second order terms in the density expansion into account, Eqs. (p3D 
and (^^, one obtains the following form of the correlation functions in the dilute // gas 
model 

D{x) = 8(n4 + -n,) I{x/p) + 9n'j^ ij}±±^^ ^(2)^^/^) 

(45) 

where I{x/p) denotes one of the possible parametrizations of the first order contribution (for 
physical reasons preferably Iph{x/p)). We stress again that in agreement with the lemma 
above the violation of (anti-) self duality, not only in instanton-anti-instanton but also in 
instanton-instanton and anti-instanton-anti-instanton superpositions, leads to the contribu- 
tions oc (?T,4 + 72.4)^. Besides these approximative solutions, there are also exact multiinstan- 
ton or multi-anti-instanton solutions of the Euclidean field equations. They are suppressed 
by higher action, and the corresponding distribution of their collective coordinates is even 
less certain within the real vacuum. Fortunately, they cannot give a contribution to Di 
due to their strict (anti-) self duality. A comparison between the leading and second order 
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(in the density) contributions to D leads to an upper bound for the packing fraction 
n\ <^ 3^r^ — TV — ■ ^ong as this holds, the density expansion of the dilute gas 
approximation should be reliable. 

With our numerical solution for jj^-* one can see from (^5|) that the instanton contribution 
to Di{x) is negative for all x, in contrast to the non-perturbative contribution as extracted 
from best fits to lattice data in This fact rises several questions. First of all, 



one might wonder whether the existing lattice data are incompatible with a negative non- 
perturbative contribution to Di, and hence are in contradiction with the instanton-gas 
model. Our experience with the existing lattice data tells, that the fits are not conclusive 
concerning the sign of the pure exponential contribution in Di. Other reasonable fits are 
possible as well. One such solution found by an own fit to the lattice data of Ref. ^ is 
shown in Fig. 2. The corresponding parameters in are 

A B 

-x~2.52xl0^, -X ^ -5.14x10^ , a ^ 0.415 , 6^ 0.307, 

^4 Ai 

AA~0.22fm, AB^0.24fm, A^ = Ai ^ oo . (46) 

This fit achieves a reasonable value /N^.f. — 3.45. 

A second question directly addresses the phenomenological consequences of negative Di. 
It was shown in |Q that D and Di define the scalar and spin-dependent potentials of heavy 
quarkonia. In functions D and Di of Gaussian shape have been used to predict the spin- 
orbit and hyperfine splitting of charmonium and bottomonium, and it was shown that only 
negative Di can reproduce the experimental situation. Therefore, it is quite possible that 
D and Di, as obtained from the instanton gas model, can generate a phenomenologically 
reasonable spectrum of heavy quarkonia. Work in this direction is planned for the future. 

The resulting combined fit of both D and Di with our model dependent input functions 
(using Iph) is shown in Fig. 3 together with the fitted exponential contributions to the lattice 
data for quenched QCD. One can see that it is still possible to imitate lattice data with D 
and Di constrained by the instanton model. The result of the best fit including the second 
order correction is the following: 
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quenched QCD 


Ir 


Iph 






0.78 
4.76 


1.29 
3.42 


1.85 
2.97 



One should note that the shape of /, I^'^\ I^'^ as functions of x/ p depends on the path 
taken for the Schwinger phase factor and on the profile of individual instantons. Concerning 
the choice of the path the straight line path (corresponding to Iph) is the most natural one 
in view of the way the correlators are measured on the lattice. As far as the profile is 
concerned, in the present paper we have chosen the exact classical solution (|lD, while its 
shape is expected to be changed due to the interactions between the (anti-)instanton and 
the vacuum (medium) on the classical and quantum level. Therefore, we can choose a more 
general instanton potential 



X 



(47) 



as proposed in with a profile function / behaving for large x like ~ exp(— mx), with 



m 



27 ^2 



TT [n^ + rii) p . This will certainly allow to describe the lattice results in a better 



way. The numerical evaluation of this effect is now in progress. 



V. CONCLUSIONS AND DISCUSSION 

We have considered the semiclassical approximation for the non-Abelian field correlators. 
The dilute instanton gas model was used. Generically, the correlators for an individual 
background configuration had to be evaluated numerically, partly in order to take into 
account the Schwinger line correctly. 

Let us briefly summarize our main results. The comparison of the considered correlators 
with lattice data at zero temperature ||^, [|1^], shows the following. The nonperturbative 
part of the correlator D can be reasonably described by the mixed I — I gas (see Fig. 3). It 
was shown analytically that selfdual configurations contribute only to the correlators D and 
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D2 in the leading order of the density expansion, while Di appears only in the second order. 
This may explain qualitatively the large ratio D/Di found in lattice simulations | 



Another interesting result is that non-positivity of Di is characteristic for the instanton 
gas model. We have demonstrated that this feature is compatible with lattice data. More- 
over, we have argued that this might be even favourable from a phenomenological point of 
view. 

The fitted values of the instanton densities turn out to be large in the case of quenched 
simulations (several instantons per 1 fm^). To describe this case consistently within the 
instanton density expansion a medium correction for the instanton profile [|l^] must be 
taken into account. In the case of full QCD (with dynamical fermions) the situation is more 
safe for the naive instanton superposition. For this case, the estimated values for instanton 
density and radius are found in the right ballpark known from instanton phenomenology, 
as well as the estimated gluon and quark condensate values extracted from the lattice data 
were consistent with QCD sumrules. 

In this paper the effect of the inclusion of parallel transporters (Schwinger lines S) into 
a semiclassical calculation (and the path dependence) has been studied for the first time. It 
is displayed in Fig. 1 that the effect (compared with a choice of path with S = 1) is of the 
order of 20 % in the region of physical interest. 

We did not discuss here the question whether instantons, which yield phenomenologically 



realistic hadron correlators [jT8[ give a complete description of the QCD vacuum configura- 
tions. The point is that the dilute instanton gas model (also in the way we have treated 
it here) does not confine. Within the vacuum correlator model, the nonzero contribution 
to the string tension which would be obtained from our correlator D, would be cancelled 
by higher correlators (see for discussion and more references). Therefore to get realistic 
vacuum there are at least two possibilities: to modify the instanton model in a way that in- 
stantons give only contributions of Gaussian nature (some Gaussian instanton ensemble) or 
to assume additional contributions, like dyons and antidyons. Work in this second direction 
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is now in progress. 
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FIGURE CAPTIONS: 



FIG. 1: The one-instanton contributions Ir,ph,s to the correlator D, according to Eq. (p6|), 
for different paths in the Schwinger hne factor S, Eq. (0). The upper sohd curve refers to 
5 = 1 for the regular gauge (anti-)instanton ( J^), the dashed line to the straight-line path 
( Iph), and the lower solid line to (anti-)instantons in the singular gauge with the infinite 
arc path as explained in the text ( Ig). 

FIG. 2: Fits to D± (upper curves) and D\\ (lower curves), as defined in (^, given in units of 
(lattice spacing- Ai,)^ versus physical distance (in units of fermi). The lattice data are taken 
from Ref. [^]. Solid lines represent the fit according to the parametrization (^,^). Dashed 
lines show our own fit (l46| ) with a negative Di, as discussed in the text. 
FIG. 3: The correlator D with the second order corrections (|45|) taken into account (upper 
solid line) and the correlator \Di\ (lower solid line) as given by the instanton gas model 
compared with the non-perturbative parts extracted from lattice data for quenched QCD p| 
(fit (|46|), upper and lower dashed line, respectively). For comparison, the long-dashed line 
shows the first order density contribution to D with the straight Schwinger line expression 
Iph- The evaluated curves correspond to the parameters n\ = 3.42 fm~^, p = 1.29Xa- 
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FIG. 1. 
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FIG. 2. 
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FIG. 3. 
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